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ABSTRACT 

In order to examine if the next generation space VLBI, such as VSOP-2 (VLBI Space Observatory 
Programmc-2), will make it possible to obtain direct images of the accretion flow around the M87 
black hole, we calculate the expected observed images by the relativistic ray-tracing simulations under 
the considerations of possible observational errors. We consider various cases of electron temperature 
profiles, as well as a variety of the distance, mass, and spin of the M87 black hole. We find it feasible 
to detect an asymmetric intensity profile around the black hole caused by rapid disk rotation, as 
long as the electron temperature does not steeply rises towards the black hole, as was predicted by 
the accretion disk theory and the three dimensional magnetohydrodynamic simulations. Further, we 
can detect a deficit in the observed intensity around the black hole when the apparent size of the 
gravitational radius is larger than > 1.5// arcseconds. In the cases that the inner edge of the disk is 
located at the radius of the innermost stable circular orbit (ISCO), moreover, even the black hole spin 
will be measured. We also estimate the required signal-to-noise ratio IZsn for achieving the scientific 
goals mentioned above, finding that it should be at least 10 at 22 GHz. To conclude, direct mapping 
observations by the next generation space VLBI will provide us a unique opportunity to provide the 
best evidence for the presence of a black hole and to test the accretion disk theory. 
Subject headings: accretion, accretion disks — black hole physics — galaxies: nuclei — Galaxyxentre — 
techniques :interferometric — relativity 
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1. INTRODUCTION 

Elucidating the nature of space plasmas and their dy- 
namical behavior in a strongly curved space-time re- 
mains one of the greatest challenges of the observational 
astronomy and astrophysics in this century. This re- 
quires a superb angular resolution, down to tens of micro- 
arcsec and, hence, has been impossible up to now with 
any instruments at any wavelengths. Surprisingly, how- 
ever, the situations is expected to be enormously im- 
proved in the near future thanks to the rapid progress of 
the very long baseline interferometer (VLBI) technique 
from the radio to sub-millimeter band. In many past 
studies, the observational feasibility of the direct imag- 
ing of the black hole shadow (or the black hole silhou- 
ette) with sub-millimeter interferometers on the ground 
has been inve stigated for the Galactic Center black hole 
Sgr A* (e.g . IFalcke. Melia fc Ago J l2000t iMivoshi et all 
2001 120071: iFish et all 120091: IDoeleman et all 120091: 
Huang. Takahashi fe Shenl I2009D and the black h ole in 
the elliptical galaxy M87 (jBroderick fe Loebll2009l) . Fol- 
lowing the same line, in this paper we elucidate the feasi- 
bility study of obtaining the direct imaging of the black 
hole with planned space- VLBI satellite. 

It is long believed that radio emission arises from 
bi-directional jets emanating from the accretion flow, 
and not from the accretion flow itself. This may not be 
the case, however, if there exist high-energy electrons 



with energy above 100 keV and if the magnetic field 
strengths are on average greater than about one percent 
of the equi-partition value. If these two conditions 
are satisfied in accretion flow, we expect significant 
radio emission not only fr om the jets but also from 
the fl ow itself (see , e.g ., iNaravan. Yi fe Mahadevanl 
119951: INaravan fe Y] I1995T ). It is true that the radio 
emission from the flow is usually overwhelmed by that 
from jets owing to much higher antenna temperature 
of the jets, but it is nevertheless possible, in principle, 
to separately detect radio emission from the accretion 
flow if observed with instrumentations with good spatial 
resolution. Future direct imaging observations of the in- 
nermost region of AGNs with VLBI will provide unique 
and excellent opportunities to m ap the accretion flow 
plung i ng into black holes (e .g.. IFalcke. Melia fc Agoll 
' 2000t IMivoshi et a ll I2004t i Takahashil 12004 12005 : 



Broderick fc LoebT 
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l200fi 
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Huan g et al 



Broderick fc Loebl l2009T ~ lNagakura fc Ta kahashi] [201 
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Such observations will also prove th e Kerr geometry 
(jBardeenl 119721: iTakahashil 12004 I2005T) predicted by the 
gravity theories, including general relativity (|Psaltisl 
I2008allbfl . 

One of the leading, on-going projects of the radio 
VLBI is the VSOP-2 with the ASTRO-G satellite, 
space VLBI program with the highest angular res- 
olution of 38 //as at the observatio n frequency of 
43 GHz (e.g. iHirabavashi et al.l 120051 iTsuboi I 120081: 
iHagiwara et al.1 120091 ) . ' The minimum antenna tempera- 
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ture which can be detected by the VSOP-2 / ASTRO-G 
is about a few times 10 8 K, which is far exceeded 
by the expected te mperature of optically th i n ac- 
cretion flow (e.g. [NaravanT Yi fc Mahadevanl H995I: 
Naravan fc Yl 119951: Esin et all 119961: iNakamura et al I 



1991 



Yuan et al 



Manmoto. Mineshige 
120031 ). 



Kusunosd 119971: 



For details, see the mission 
web site (http://www.vsop.isas.ac.jp/vsop2e/). 
There is another on-going space- VLBI mis- 
sion called RadioAstron (for details, see, 
http://www.asc.rssi.ru/radioastron/). We have 
mainly used the parameters of the space- VLBI satellite 
in VS0P-2/ASTR0-G as one example and the results 
presented in this paper can be applied to other similar 
space- VLBI missions. 

The existence of a black hole in an accretion flow 
or a disk will manifest itself as a region with a deficit 
intensity, sometimes called as "black hole shadow" or 
"black hole silhouette", in the observed intensity map. 
The properties of the black hole shad ow were inves- 
tigated for optically thi ck disks (e.g. iLuminetl 119791 : 



Fukue fe Yokovamal 19881) a nd for optically thin flow (e.g. 



Falcke. Melia fe Agolll2000D . The angular size and precise 



shape of the black hole shadow depend on the distance, 
black hole mass, space-time geometry, and the optical 
dept h of the accret i on flow at the observed frequencies 
(e.g. iBardeenl 119721 ILuminetl [l97l iFukue fc Yokovamal 
Il988t lTakahashill2004 120051) . 

There are good targets of supermassive black holes 
suitable for the VLBI imaging observations. For the ob- 
servations at 43 GHz, the best target will be M87, which 
has the second largest angular size of the gravitational 
radius, r g /D ~ 1-3 /ias. Here, r g = GM/c 2 , and M 
and D are the mass and the distance of the black hole, 
respectively. According to the relativistic ray-tracing cal- 
culations under the Kerr geometry, the maximum size of 
the shadow ran ges between 4-15 r g /D (sec, Fig. 2 in 
lTakahash3l2004D which corresponds to 6-24 /zas for M87 
when il/ = 3x 10 9 M Q and D = 16 Mpc are assumed 
where M Q is the solar mass. Note that the black hole in 
the Galactic Center (Sgr A*), the largest black hole in 
terms of the apparent angular size oir g /D ~ 5 /xas, is not 
suitable for the observations at 43 Gflz, since the image 
of the accretion flow around Sgr A* will be totally washed 
out at low frequencies, < 100 ~ 150GHz, because of in- 
terstellar scattering. In fact, a rather broadened radio 
morphology with an elliptic shape was already obtained 
at the observed frequency 86 GHz (e.g., Shen et al. 2005, 
see also Doeleman et al. 2008 for the observations at 230 
GHz). 

In this study, we investigate how the accretion flow 
around the black hole in M87 will be observed and 
mapped by technically feasible space- VLBI observations 
at 43 GHz and 22G Hz. For the bolomet r ic luminosity of 
~ 10 41 erg s" 1 dBiretta et al.l ll991t iRevn olds et al.l 
119961: iDi Matteo et al.l 12003ft and the black hole mass 
of M ~ 3 x 10 9 Mo, the Eddington ratio of M87 is 
about L/Le ~ 10~ 6 , much less than unity, indicat- 
ing that the accretion flow i s likely to be radiatively 
ineffi cient (see Chapter 9 of iKato. Fukue fc Mineshigg 
120081 ) . The observed multi- wavelength continuum spec- 
trum of M87 from radio to X-rays can be explained 
in terms of the RIAF (radiatively inefficient accre- 



tion flow) model (or its varian t, advection-dom i nated 
accretion flow, ADAF model) (|Di Matteo et al.l 
Wang et al.ll20 08: Bro derick fc Loebll2009t l~ et al.l 



2003 



2009; 



Nagakura fc Ta kahashi 20101) . Further, the powerful, rel- 
ativistic jets are observed in M87, which seem to be 
produced in the innermost part of the accretion flow. 
To understand the rapid TeV 7-ray variability in M87 
jets detected by the High Energy Steroscopic System 
(HESS, Aharonian et al. 2006), it is sugges ted that the 
black hole in M87 may be rapidly rotating (jWang et al.l 
I2008D . On this ground, we calculate the expected im- 
ages of the central region of M87 by the radiative trans- 
fer calculations in the Kerr space-tim e, basically adopt - 
ing the same method as tha t used in iTakahashil (120041) . 
ITakahashi fc Watarail (|2007t) and ITakahashi fc Haradal 
( 20101) except that we employ the RIAF model, instead of 
the standard disk model, in the present study. The con- 
tinuum emission at around 43 GHz will be of Rayleigh- 
Jeans spectrum because of the self- absorption of the syn- 
chrotron radiation by the thermal electrons. We, there- 
fore, assume that the accretion flow is optically thick to 
absorption at 43 GHz. 

The plan of this paper is as follows: In §2 we explain 
the methods and models in this work. The results of 
numerical simulations are given in §3. The important 
issues are discussed in §4 and the conclusions are given 
in §5. 

2. METHODS AND MODELS 

We first need to specify the black hole mass M and the 
distance D to M87, but there are uncertainties in the ob- 
served values of M and D. In FigfTJ we summarize the 
past reports for the mass and the distance measurements. 
In this figure, we plot three data for the mass measure- 
ments and four data for the distance measurements with 
uncertainties indicated by the shaded region (for the de- 
tails of the data, see the caption of Figjj). In this study, 
we adopt three datasets of combination of M and D de- 
noted as A, B and C in FigfTJand these values are sum- 
marized in Table [T] In contrast, the black hole spin is not 
so well constrained by the observations. Therefore, we 
consider two extreme cases: a black hole with spin pa- 
rameter a = (i.e., a Schwarzschild black hole) and one 
with a — 0.998 (i.e., an extreme Kerr black hole) where 
a = J I GM 2 c~ x where J is the angular momentum of the 
black hole. 

Next, we need to prescribe electron temperature dis- 
tribution. According to the RIAF model, ion temper- 
atures are very high, close to the virial, and are sev- 
eral times 10 12 (r/r g )~ 1 K near the black hole where 
r is the radial distance from the center, whereas the 
electron temperature will be much less and saturated 
around a few times 10 9 K (which is about m e c 2 /k with 
m e and k being the electron mass and the Boltzmann 
constant, respectively) due to the enhanced coolin g in 
the relativistic regimes (e.g., INaravan fc YH II995L see 
Kato, Fukue, & Mineshige 2008 for a review and ref- 
erences therein). Note that the accretion disk corona 
model also gives similar electr on temperature profiles but 
with slightly lower values (e.g,. ILiu. Mineshige fc Ohsugal 
I2003D . As a result, the electron temperature is spatially 
uniform close to the black hole while the ion temper- 
ature increases in the central region of the disk. Such 
features are confirmed by three-dimensional MHD sim- 
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TABLE 1 

The black hole mass of and the distance to M87 adopted in this 
study. The corresponding apparent angular sizes are also 
listed. 



model 


A 


B 


c 


black hole mass M [10 y M Q ] 


6.4 


3.0 


1.7 


distance to M87 D [Mpc] 


17.9 


16.0 


17.9 


angular size corresponding to lOGM/(c 2 D) [/lias] 


35.0 


18.5 


9.3 
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Fig. 1. — Diagram of the black hole mass, M, and the dis- 
tance, D, with the observational data points of M87. Lines of 
WGM/(c?D) = 10/^as, 20/jas, 30/xas, 40/xas and 50/ias are also 
indicated. The ranges of the black hole mass of and the distance 
to M87 reported in the past studies are plotted by the shaded re- 
gion. For t he bla ck hole m ass, we use M = (2. 4 ± 0.7) x 10 9 M o 
[Ford et al.l fl99l)(F+94); IHarms et alJ l(l994|)(H+94)], (3.2 ± 
0.9) X 1 O 9 M jMacchetto et all pTj97T l (M+97)] and (6.4 ± 0.5) X 
10 9 Mq Gcbhardt & Thomas (2009) (GT09)]. The mass in GT09 
is calculated by assuming D = 17.9 Mpc. For t he distance to 
M87, we use D = 14.7±1.0M pc rjacobv et alJ [119901) (J +90)], 16.75 
Mpc I Whitmore et al.l lli995l) (W+95)1 , 16.0±1.9 Mpc IMacri et al.l 
l |1999M M+00)] and 17.0 ± 0.3 Mpc jTonrv et all H200U )(T+01)j. 
The values of M and D adopted in the present study are denoted 
by A, B and C, respectively (see also Table [TTl . 



ulation (e.g.. lOhsuga. Kato fc Mineshigell2005[) . The ex- 
pected electron temperature of > 10 9 K in the cen- 
tral part of the flow certainly exceeds the detection 
limit of VSOP-2 / ASTRO-G at 22 GHz and 43 GHz. 
Hence, we prescribe electron temperature profile as 
T c = 5 x 10 9 (r/Vb) -p K where p and r& are param- 
eters. In this study, we calculate three cases of p: 
p = 0.1, 0.5 and 1.0. Note that the last model 
has the same temperature profile as that of ions. As 
for the parameter rb, we determine the value so as 
to reproduce the observed energy spectrum, finding 
rb ^ 100r g . Similar electron t emperature profile s were 
used in the past study (e.g . iBroderick fc Loebl 120091 : 
INagakura fc Takahashil I2010D and successfully repro- 
duced the observed energy spectrum and images. At 22 
GHz and 43 GHz, the energy spectrum is well explained 
by t hermal synchrotron emissio n (e.g. iDi Matteo et al ' 



20i 



2003HWang et a l. 2008; iLi et alj|2009t IBroderick k Loeb 



120091: INagakura fc Takaliaihll2^Tol) . 

When the radiative cooling is not so efficient above 
~ 10 10 K and/or the energy transportation from ions 
to electrons is very efficient, electron temperatures can 
increase with decreasing r. So, the radial profile of the 
electron temperature contains the important information 
about the electron heating mechanism of the accretion 
flow. In this study, we show that the observed images 
obtained by the space- VLBI observations can put useful 
constraints on the electron temperature profile. 

Further, we need to prescribe the inner edge and the 
velocity pattern of the accretion disk. We consider two 
cases for the inner edge of the luminous part of the ac- 
cretion disk: one is at the innermost stable circular orbit 
(ISCO), nsco 5 and the other is at the event horizon, 
rjj- In this study, we assume the sub-Keplerian velocity 
profile with the angular velocity given as 
and the radial component of the four velocity given as 
u r = Up F /y/2 where and Up F are the Kcplcrian an- 
gular velocity and the radial four velocity component of 
the freely falling motion with zero angular momentum 
at infinity, respectively We expect that the accretion 
disk in M87 has the sub-Keplerian velocity distribution 
near the black hole. As for the cases of different veloc- 
ity profiles (such as freely falling motion and Keplarian 
rotation), see discussion in §4. The viewing angle, i, 
between the direction of the observer and the rotation 
axis of the disk is assumed to be i = 45°. Although in 
the p ast studies the inclination ang le of i = 30° is used 
(e.g., INagakura fc Takahashil 120 10D . we have confirmed 
that the main results in this study do not change even if 
i = 30° is used. Based on the assumptions described 
above, we perform the general relativistic ray-tracing 
calculations in the Kerr space time; that is, we fully 
incorporate the relativistic effects, such as the frame- 
dragging, the gravitational redshift, the bending of light 
and the Doppler boosting. Our ca lculation methods are 
prescr i bed in our past papers, i . e. iTakahashi &; Watarail 
(120071). INagakura fc Takahashil (|2010D and Appendix of 
ITakahashi fc Haradal (|2010D . 

After calculating the simulation images of the inten- 
sity /„ with observed frequency at 22 GHz and 43 GHz 
by assuming the thermal synchrotron spectrum, we cal- 
culate the expected images smeared out by the spatial 
resolution of ~ 38/i arcseconds with the realistic u-v cov- 
erage data created by the simulation tool, th e Astronom- 
ical Radio Interferometer Simulator (ARIS; lAsaki et al.1 
12001 . In the left panel of Fig [2j we show a sample 
of the u-v coverage of the VSOP-2 (ASTRO-G satellite 
and VLBA) for the observation of M87 within one day. 
By using this u-v coverage, the synthesized dirty beam 
is calculated as shown in the right panel of Fig [2] The 

We thank Y. Asaki for providing ARIS. 
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Fig. 2. — A sample of the u-v coverage of the VSOP-2 for the observation of M87 within one day {left) and the corresponding synthesized 
dirty beam for model B [right). 



TABLE 2 

Models 



model 


M&cD 


a 


i 


rin 


rout a 


Ibh a 


Asymmetry 


Deficit 


Fig 


p = 0.5 A(a)p05 


A 





45" 


nsco 


66 


20 


Yes 


Yes 


3, 4 


A(b)p05 


A 





45° 


r H 


66 


20 


Yes 


Yes 


3, 4, 7 


A(c)p05 


A 


0.998 


15° 


r H 


66 


20 


Yes 


Yes 


3, 4 


B(a)p05 


B 





45° 


r ISCO 


110 


26 


Yes 


Yes 


3, 4 


B(b)p05 


13 





15° 


m 


110 


21 


Yes 


Yes 


3, 4, 5 


B(c)p05 


B 


0.998 


15° 




110 


25 


Yes 


Yes 


3, 4 


C(a)p05 


C 





45° 


nsco 


220 


38 


Yes 


Yes 


3, 4 


C(b)p05 


C 





45° 


■r H 


220 


39 b 


Yes 


No 


3, 4 


C(c)p05 


c 


0.998 


45° 


r H 


220 


38 b 


Yes 


No 


3, 4 


p = 0.1 A(a)p01 


A 





45° 


nsco 


66 


90 


Yes 


Yes 


4 


A(b)p01 


A 





45° 


m 


66 


90 


Yes 


Yes 


1 


A(c)p01 


A 


0.998 


45° 




66 


90 


Yes 


Yes 


1 


B(a)p01 


B 





45° 


nsco 


110 


110 


Yes 


Yes 


4 


B(b)p01 


B 





45° 


■r H 


110 


110 


Yes 


Yes 


1 


BfclpOl 


B 


0.998 


15° 


r H 


110 


110 


Yes 


Yes 


1 


C(a)p01 


C 





45° 


nsco 


220 


130 


Yes 


Yes 


1 


C(b)p01 


C 





45° 


Hi 


220 


130 


Yes 


Yes 


1 


C(c)p01 


c 


0.998 


45° 


rH 


220 


130 


Yes 


Yes 


4 


p = 1.0 A(a)pl 


A 





45° 


nsco 


66 


16 


Yes 


Yes 


4 


A(b)pl 


A 





45° 


m 


66 


15 


Yes 


Yes 


1 


A(c)pl 


A 


0.998 


45° 


r H 


66 


15 


Yes 


Yes 


1 


B(a)pl 


B 





45° 


nsco 


110 


19 


Yes 


Yes 


4 


B(b)pl 


B 





15° 




110 


18 b 


Yes 


No 


1 


B(c)pl 


B 


0.998 


15° 


r H 


110 


17 b 


Yes 


No 


1 


C(a)pl 


C 





45° 


nsco 


220 




Yes 


No 


1 


C(b)pl 


C 





45° 


m 


220 




marginal 


No 


1 


C(c)pl 


c 


0.998 


45° 


TB 


220 




marginal 


No 


4 


p = 0.5 B(a)p05il5 


B 





15" 


nsco 


110 


26 


Yes 


Yes 


6 


B(b)p05il5 


B 





15° 




110 


21 


Yes 


Yes 


6 


B(c)p05il5 


B 


0.998 


15° 


rs 


110 


25 


Yes 


Yes 


6 


a In units of r K = GM/c A . 


















b Marginal value. 
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scale of the dirty beam shown in this figure is calculated 
based on model B. 

The calculation method of the observed image smeared 
by the spatial resolution of the VSOP-2 / ASTRO-G ob- 
servations is as follows: 1) From the dirty beam, the 
primary lobes of the clean beam can be calculated by 
fitting the primary component of the dirty beam by 
the elliptical gaussian. 2) We calculate the visibility 
V(u,v), that is, the Fourier components of the elliptic 
gaussian beam S(u, v) and the theoretical images. 3) 
The map of the spatially smoothed intensity I u is cal- 
cula ted as Ijx,y) = f f e 2ni ( xu +y v *> S(u, v)V(u, v)dudv 
(e.g. iThompson et al.ll200lL hereafter, TMS). The image 
produced based on this procedure basically corresponds 
to the image produced by the CLEAN deconvolution al- 
gorithm if the algorithm works well and if the effects of 
the side lobe can be completely removed, i.e. no decon- 
volution errors (see, Sec 4 for important discussion on 
the data deconvolution method for the image containing 
the dark region around the black hole). In the CLEAN 
deconvolution algorithm (e.g., Sec 11.2 of TMS) which is 
usually used in the radio data analysis, the clean map is 
calculated by first extracting the dirty beam components 
with some factor from the dirty map and then replac- 
ing the dirty beam with a gaussian (or similar functions 
that are free from negative values). If this procedure 
works well, the effects of the side-lobe patterns can be 
removed and the resultant images made by the CLEAN 
deconvolution algorithm corresponds to the image which 
is obtained from the original image smeared with the 
gaussian. Therefore, the image produced by the calcula- 
tion method used in this study (as denoted above) corre- 
sponds to the images with no deconvolution errors and no 
effects of the sidclobe in the CLEAN deconvolution algo- 
rithm. About the possible problems in the application of 
the CLEAN deconvolution algorithm to the images with 
the dark region around the black hole, see the discussions 
in Sec 4. 

In this study, we calculate the smeared images for 30 
models (for model parameters see Table [2]). We first 
prescribe the outer radius of (the luminous part of) the 
accretion disk as follows. According to the past observed 
images of M87 by the radio interferometers the size of 
an emitting region is a bout a few to several hundreds 
of mi cro arcseconds fe.g. Uunor et al.lll999tlDodson et al.l 
I2006T ) . Note, however, that the actual length scale should 
depend on the distance, D, for a fixed value of the an- 
gular size, and the gravitational radius depends on the 
mass, M. In order that the emitting region size should 
be consistent with the past observational data, we adopt 
r out = 66, 110, and 220 (r g ) for Models A, B, and C, re- 
spectively (see Table [2]). As to the inner radius and the 
spin parameter, we consider three sets of their combina- 
tions: (a) (a, r in )= (0, risco), (b) (0, r H ), and (c) (0.998, 
?"h)- From the model names, we can get the information 
of the assumed parameters, which include M, D (from 
A, B or C, sec Tabled]), a, In [from (a), (b) or (c)], and 
p (from pOl, p05, pi). 

3. APPEARANCE OF THE ACCRETION DISK AROUND A 
BLACK HOLE 

3.1. Overview 

Fig |3] shows the results of the observed intensity distri- 
butions for the case of p = 0.5 (i.e., T c cx r -0 5 ). From 



the left to right panels, wc show the theoretical (un- 
smoothed) images (first column), images smeared with 
the u-v coverage of the space- VLBI observations based 
on models A, B and C in Table [1] (second, third and 
forth columns), and the images smeared with an ellipti- 
cal be am size of 0.3 mas x 0.2 mas as in the past observa- 
tions ([Dodson et al.|[2006() (fifth column) . From the top 
to bottom panels, we show the results for the above cases 
(a) (top), (b) (middle) and (c) (bottom). In this figure, 
the normalized intensity variation along the horizontal 
line passing the center is also shown for each image ( blue 
solid line). For other informations, see the caption of 
FigH 

For all the cases of the space- VLBI observations (sec- 
ond, third and forth columns), the asymmetry of the ob- 
served intensity with respect to the vertical line pass- 
ing the intensity peak can be clearly seen, while in the 
smeared images by the resolution of the past observa- 
tion (fifth column) such asymmetry can not be seen. 
In all the cases of the space- VLBI observations except 
models of C(b)p05 and C(c)p05, furthermore, an under- 
luminous region (i.e. deficit of the observed intensity) 
around the central black hole can be seen; that is, sec- 
ond peak of the intensity can be clearly seen. The size 
of the under-luminous region is nearly the same as that 
of the beam size (gray ellipse) for models of A(a)p05, 
A(b)p05, A(c)p05 and B(a)p05 and is smaller than the 
beam size for models B(b)p05, B(c)p05 and C(a)p05. On 
the other hand, the intensity contrast amounts to ~ 40% 
for A(a)p05, - 30% for A(b)p05, A(c)p05, and B(a)p05, 
~ 20% for models B(b)p05 and B(c)p05, and - 10% for 
model C(a)p05. In terms of the black hole spin, we can 
not see the difference between a = and 0.998, in the 
smeared images of the cases of n n = rn (compare middle 
and bottom), while we can see the difference with ~ 10% 
difference of the intensity in the case of rj n = risco 
(compare top and bottom). Here, note m ~ nsco for 
a = 0.998. 

Here, we define the size £bh as the length between 
the two peaks of the intensity profile in units of GM/c 2 , 
which gives the size of the under-luminous region around 
the central black hole, i.e. from the value of £bh, the 
spatial size containing the black hole can be measured. 
When there is only one peak in the intensity profile, the 
size £bh can not be defined. In Table [2] we summa- 
rize the values of £bh for all models. In the case of Sgr 
A*, the VLBI observations put constraints on the size of 
the region containing the black hole as ~ 25.2 GM/c 2 
from the size me asurements of the observed intensity 
(|Shen et al J 120051 ). In Table [H we can see that mod- 
els A(a)p05, A(b)p05, A(c)p05, B(a)p05, B(b)p05 and 
B(c)p05 give constraints of £bh around £bh ~ 25 GM/c 2 
corresponding the size obtained in Sgr A*. In all mod- 
els with p = 0.1, wc find £bh ^ 100. In models of 
A(a)pl, A(b)pl, A(c)pl, B(a)pl, B(b)pl and B(c)pl, on 
the other hand, the size £bh is smaller than < 19 GM/c 2 
which gives a more stringent constraint than the case of 
Sgr A*. If this is the case, the space- VLBI observations 
can put useful constraints on the size £bh of the region 
containing the black hole. In Table [2j we also indicate 
whether the asymmetry (ninth column) and the deficit 
(tenth column) in the intensity profile can be seen or not 
for all the calculated models. Although the asymmet- 



6 



Takahashi and Mineshige 



original 




100GM/{c 2 d) 350/ias 



Fig. 3. — Images of an accretion flow surrounding a black hole with the characteristics of M87. The emitting accretion flow is assumed 
to have a sub-Keplarian velocity field with an electron temperature profile of T e = 5 X 10 9 K (r/r^)~ p with p = 0.5 and = fOO r g . From 
the left to right panels, the theoretical (unsmeared) images (first column), the images smeared with the u-v coverage of the next-generation 
space- VLBI observations [with the beam size of 38/^asx 120/xas (gray ellise)] for models A, B and C ( in T able U (s econd, third and forth 
columns), and the images smeared by the beam size of the past observation (with 0.3 masx0.2 mas, UDodson et al.ll2006ft (fifth column), 
respectively. The black hole is either (a, b) non-rotating (a = 0) or (c) maximally rotating (a = 0.998). The viewing angle is fixed to be 
i = 45° in this figure. The inner edge of the accretion flow is either at (a) the radius of the innermost stable circular orbit (ISCO) nsco ° r 
(b, c) at the event horizon rjj. The left panels show the results based on models A(a)p05, A(b)p05, and A(c)p05, while the right panels arc 
based on models B(a)p05, B(b)p05, and B(c)p05. These model parameters are summarized in Table[2] The rotation axis of the accretion 
flow and the black hole is also shown (green dotted lines) in the left panels. The normalized intensity variation along the horizontal line 
passing the black hole center are overlayed (blue lines). The vertical and horizontal width of images is 100 GM/(c 2 D) which corresponds 
to about 295/xas (A), 185^tas (B), and 93/^as (C). 



ric signature can be seen (or marginally seen) for all the 
models, the deficit can not be seen for some models. 



3.2. Temperature profiles 

Next, Fig 2] compares the cases of different electron 
temperature profiles: p — 0.1 (green solid line), 0.5 (blue 
solid line) and 1 (pink solid line). In this figure, we show 
the normalized intensity variation along the line passing 
the center for 27 models listed in Tablc[5] All the lines for 
p = 0.5 are the same as those of Fig[3]and the lines for p = 
0.1 and 1 are calculated in the same way as those in Fig[3J 
For all the cases of the space- VLBI observations (second, 
third and forth columns) except for models of C(b)pl and 
C(c)pl, the asymmetry of the observed intensity with 
respect to the vertical line passing the intensity peak can 
be clearly seen, while in the cases of C(b)pl and C(c)pl, 
only the asymmetry can be marginally seen. 

On the other hand, we can clearly see the under- 
luminous region around the black hole for all models ex- 
cept C(b)p05, C(c)p05, B(b)pl, B(c)pl, C(a)pl, C(b)pl 



and C(c)pl. In general, the steeper the radial tempera- 
ture profile is (e.g. p = 1), the more difficult it becomes 
to see the asymmetry of the intensity and the under- 
luminous region around the black hole. This is because 
the size of the effectively emitting region for p = 1 is 
smaller, compared with other cases (p = 0.1 and 0.5) and 
so the luminous part of the image tends to be smeared 
out. The observational feasibility of the effects of the 
black hole spin is basically the same as that in the case 
of p = 0.5 described above, i.e. in the smeared images 
of the cases with r- ln = rjj (compare middle and bottom), 
we cannot see a difference between the cases with a = 
and 0.998, while in the case with r m — nsco (compare 
top and bottom) we can see. 

3.3. Observational errors 

So far, we implicitly neglected the effects of the noise 
introduced by the observations. Even when the wide 
u-v coverage is achieved by the space- VLBI satellite, 
the observational features investigated above if noise is 
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Fig. 4. — The normalized intensity variations along the line passing the black hole center for p = 0.1 (green lines), 0.5 (blue lines), and 
1 (pink lines), respectively. All the lines for p = 0.5 are the same as those in Fig|3] The normalized intensity variations are shown for the 
theoretical (unsmeared) images (first column), the images smeared with the u-v coverage of the next-generation space- VLBI observations 
for models A, B and C (in Tablc[l} (second, third and forth columns), and the images smeared with the resolution of the past observation 
resolution, respectively The lines of p = 0.5 are the same as those in Fig. \3\(fifth column). The vertical width is 100 GM/(c 2 D) which 
corresponds to about 295/^as (A), 185/^as (B), and 93^tas (C). See Table [5] for model parameters of the plotted 27 models. 



large, will not be completely obtained. The noise is 
caused by several factors; e.g. deviations of the antenna 
surface from the ideal profile (for details, see, TMS). 
The complex visibility V is deviated by the noise as 
Z = V + e where Z is the measured visibility (includ- 
ing the effects of the noise) and e is the noise compo- 
nents (see, Fig 6.8 in TMS). The effects of the noise 
can be included by assuming Gaussian noise of stan- 
dard deviation <7 no isG, and the probability distributions 
of the measured visibility is given in Sees 6.2 and 9.3 
of TMS. In the space- VLBI observations, the standard 
deviation cr no ise is calculated as the flux density calcu- 
lated by using the SEFD Se of the space- VLBI satellite 
(S E pacc ) and the ground-based interferometer (S'| round ) 

as o-noiso = [5E Pacc S'| round /(2A^T a )] 1 / 2 /?7Q where Av is 
the bandwidth, r a is the data averaging time and tjq is 
the efficiency factor (see, e.g. Sec. 6.2 of TMS). In this 
paper, we use the signal-to-noise (SN) ratio IZsn defined 
as 7?. S n = Oscurcc/o-noisc where a- sourC c is the spatial aver- 
age value of the source flux density in the region within 
the outer disk radius r out . 

Based on the probability distribution of the noise (in 
Sec 9.3 of TMS) and the assumed value of TZsn, the ef- 
fects of the noise can be included. In this paper, we 
take the following procedure to include the noise effects: 



1. First, we calculate the average flux density <7 S ourcc 
for each model in Table 2. 

2. Next, we prescribe a single value of IZsn, e.g. 
"^•sn = 50, and calculate the standard deviation of 
the noise <7 n oisc(= ^source /7?-sn) based on the cal- 
culated value of CT sourco and the assumed 72-sn- 

3. Here, we injected a single noise value and input the 
noise effects into (u, w)-points. For the single value 
of (Tnoiso calculated in the last stage and the gaus- 
sian probability distribution function of the noise 
(given in Sees 6.2 and 9.3 of TMS), we can pro- 
duce random gaussian noises. In this stage, we 
randomly input the gaussian noise on each (u, v )- 
point of the measured visibility which is calculated 
from the theoretical image by the Fourier transfor- 
mation, i.e. we have the noisy visibility map in this 
stage. It is noted that at this stage since the noise 
is randomly produced for each visibility point, the 
noise values arc different at different (u, v) -points. 

4. From the noisy visibility calculated at stage 3, we 
convolve the sampling function calculated from the 
u-v coverage of the space- VLBI satellite and the 
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ground-based interferometers. We then calculate 
noisy image by the inverse Fourier transformation. 

5. We make ten independent noisy images in this 
stage, i.e. stages 3 and 4 are repeated for ten times 
based on the different gaussian random noise pat- 
terns. 

6. Finally, we average the ten noisy images. 

In this procedure, we calculate the noisy image from the 
theoretical original image for an assumed value of TZsn- 
It is noted that the value of 72-sn represents the SN ratio 
for the averaged value of the source flux density (from 
the definition of 72-sn) and the actual SN ratios at differ- 
ent points of the image (or the visibility) are different. 
In general, around the most luminous part of the image 
where the value of the flux is larger then the averaged 
flux, the SN ratios are much higher than 7?-sn which rep- 
resents the average SN ratio. 

In the left panel of Fig[SJ we show the normalized inten- 
sity variations for model B(b)p05 with noise of 7?-sn = 10, 
and 50 and without noise. These curves with noise are 
calculated by the average of the ten curves with different 
noise patterns. In this way, we find that a deficit of the 
intensity can be measured for 7?-sn ^ 10. In the curves of 
= 10 and 50, another deficit of the intensity caused 
by the noise can be seen around x = 40. The depths of 
these deficits (around x ~ 40) are comparable to those 
of the deficit in the center caused by the existence of the 
black hole. We also calculated the case for 7\Lsn = 100 
and confirmed that the deficits around x = 40 disappear. 

Here, for the ASTRO-G satellite, we estimate the value 
of 5g pacc satisfying the condition of 7£sn > 10 or > 100. 

By using the values of 

f source — 

1 Jy for M87, Sf ound = 
1500 Jy for the Very Large Baseline Array (VLB A), tjq = 
0.7, Ai/ = 256 MHz and r Q = 60 s, the SEFD required 
for the observations satisfying 7?-sn ^ 10 is 5g at < 
1.0 x 10 5 Jy, and for ft SN > 100, the required SEFD is 
^satellite < i.q x 10 3 Jy. At 43 GHz, the nominal and 
the possible worst values of g| atelllte are 28000 Jy and 
190000 Jy for the ASTRO-G satellite. Q Therefore, we 
can expected the SN value around 7£sn ~ 10 ~ 100 for 

(~ 28000 Jy). For the worst 
(~ 190000 Jy), the SN value is TZ SN < 7, 
and then useful images cannot be obtained. 

We also investigate the effects of the noise for the case 
of 22 GHz for model B(b)p05 (see the right panel of Fig 
[5]). The asymmetric brightness profile can be seen for 
the cases of 72-sn = 50 and 100. From these curves, 
we can marginally put a constraint on the size £bh as 
£bh ~ 30 GM/c 2 for TZ SN ~ 100 and - 35 GM/c 2 for 
IZsn ~ 50. In these cases, however, we cannot see a 
deficit of the observed intensity around the black hole. 
In the same way as we did for 43 GHz, we estimate the 
value of 5^ paco satisfying the condition of TIsn }t 50. 

By using the values of 5| round = 500 Jy for the VLBA, 
r a = 120 s and the same values for other parameters, the 
SEFD at 22 GHz required for the observations satisfying 
^sn > 50 becomes S|f tcllitc < 2.4 x 10 4 Jy. At 22 GHz, 
the nominal and the possible worst values of ^gs-telhte are 



the nominal value of S|, atelllte 
value of 5f tellite 



5000 Jy and 8200 Jy for the ASTRO-G satellite. Even 
for the worst case, we can expected the SN value larger 
than 100, i.e. TvLsn > 100. Therefore, the observations 
at 22 GHz by the VSOP-2/ ASTRO-G will be able to 
detect the asymmetric brightness profile of the accretion 
flow around the black hole. This means that such ob- 
servations will give a new evidence of the existence of a 
black hole within the size of £ B n ~ 30 GM /c 2 in M87. 

4. DISCUSSION 

In this study, we aim at extracting the essential fea- 
tures of the M87 image to be observed with a techni- 
cally feasible space- VLBI satellites. For the calculation 
purpose, we have made several assumptions and simpli- 
fications. Here, we discuss how realistic our numerical 
treatments are. 

As is well known, bright jets are emerged from the 
M87 black hole. In this study, we implicitly neglected 
their contribution to the radio emission at 43 GHz and 22 
GHz. We also neglected the nonthermal electron compo- 
nent. These treatments can be justified, since the outflow 
component can be n egligible at 43 GH z , compared with 
the d i sk component (JJunqr et al.lll999t iDi Matteo et al.l 
120031: iDodson et all [2006]), although it should be in- 
cluded in the othe r energy bands (jBroderick fc LoeO 
l2009tlLi et all[200l . Further, the observed energy spec- 
trum can be fitted by the thermal synchrotron as was 
already mentioned in the previous section. Furthermore, 
since the electromagnetic radiation from the jets is likely 
to be generate by internal shocks, which take place at a 
few tens to a few hundreds of r g apart from the footpoint 
of the jet; that is, the luminous parts of the jets should be 
very distant from the black hole. It will be thus easy to 
discriminate the radio emission from the disk from that 
from the jets. 

In this study, we only consider the sub-Keplcrian veloc- 
ity distribution since the flow is very likely to be RIAF, 
which typically shows sub-Keplerian ro tation. In the 
past studies (e.g. iBroderick fc Lo eb 2009) , the Keplerian 
rotation model is considered. Since the Doppler boost- 
ing effects become much stronger in the case of the Ke- 
plarian rotation, we expect that the asymmetry of the 
intensity variation shown in Figs [3] and |4] will be more 
prominent. In the case of the freely falling motion, con- 
versely, the asymmetry of the intensity variation will be 
less pronounced, compared with the sub-Keplerian cases 
considered in this study. 

In this study, we also neglected the effects of a finite 
disk thickness. Since in the case of M87 the viewing angle 
is not large, i.e. i < 45°, the effect of the disk thickness is 
negligible and we have confirmed this by including the ef- 
fects of the thickness for some samples. Concerning the 
viewing angle, some past st udies used more smaller v alue, 
such as i = 15° ~ 20° (e.g. IBroderick fc Loebl[200l and 
references therein) . We have confirmed that the main re- 
sults in this study will not be changed even for such small 
viewing angles. In Fig[6j we plot the normalized intensity 
variations along the line passing the black hole center for 
i = 15° [models B(a)p05il5, B(a)p05il5 and B(c)p05il5] 
and 45° [models B(a)p05, B(a)p05 and B(c)p05]. The 
curves for the cases of i — 45° are same as those plot- 
ted in Figs |3] and |4] From these plots, we can clearly 
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Fig. 5. — Normalized intensity variations along the line passing the black hole center for model B(b)p05 without noise (the blue lines) 
and with noise (the brown lines) of the signal-to-noise ratio = 10 and 50 for the observed frequency 43 GHz (left panel) and 7?-sn = 50 
and 100 for 22 GHz (right panel). For the case of 43GHz, the curve of no noise is same as the curves of the model B(b)p05 shown in Figs 
[3]andH 
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Fig. 6. — Normalized intensity variations along the line passing the black hole center for i = 15° [models B(a)p05il5, B(a)p05il5 and 
B(c)p05il5] and 45° [models B(a)p05, B(a)p05 and B(c)p05]. Model parameters are given in Table[2] 



see asymmetric signatures and central dark regions for 
the cases with i = 15° as in the cases of i = 45°. For 
the cases with i ~ 15°, the size £bh which is the length 
between the two peaks of the intensity profile is also cal- 
culated and is summarized in Table [5J The differences 
in the intensities of the two peaks in the cases of i = 15° 
are smaller then those in the cases of i = 45°. This is be- 
cause in the cases of i — 15° the Doppler boosting effects 
become smaller than the cases of i — 45°. For the cases 
with i = 15° ~ 45°, we expect similar features (e.g. the 
asymmetric intensity profile, the central dark region) as 
shown in Figs [3] and [4] 

As shown in the previous section, the effects of the 
black hole spin can be seen by the space- VLBI observa- 



tions at 43 GHz for the case of rj n = ngco as far as 
the stationary accretion disk model is concerned. How- 
ever, the realistic accretion flow shows significant time 
variability. For the black hole in M87, the Keplerian 
rotation timescale, T, at nsco is T ~ 9.2 x 10 4 s 
[(nscoAg) 3/2 + a] [M/ (3.0 x 10 9 M Q )] for a distant ob- 
server. This timescale corresponds to 16 d and 2.1 d 
for a non-rotating and a maximally rotating black hole, 
respectively, and are both longer than the typical dura- 
tion of a VLBI observation. Thus, the intensity inhomo- 
geneities in the M87 accretion disk and the time variable 
signatures (i.e., flares) will be detected by continuous ob- 
servation with the space- VLBI satellites. If such observa- 
tions are performed, the black hole spin will be measured 
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Fig. 7. — Normalized intensity variations along the line passing 
the black hole center for the cases using the dirty beam and the 
elliptic gaussian beam. 

from the variation timescale of the observed flux, and the 
sp ace- VLBI observati ons will be able to confirm the claim 
bv lWang etlll (120081 ) that the M87 black hole is rapidly 
rotating. Moreover, such inhomogeneity will produce a 
transient deficit of the observed flux, which is similar to 
that around the black hole. However, the former will be 
time variable while the latter will be stable both in time 
and position. It is thus easy to discriminate one from the 
other. We will investigate the details of these topics in 
near future. 

We should also discuss what information we can ex- 
tract if no deficits nor asymmetry were detected in the 
observed intensity profile. In such cases, it is expected 
that the velocity patterns should nearly be free-fall as 
long as the disk is luminous, or that only the jet compo- 
nent can be seen and the accretion disk is not luminous 
due to a low mass accretion rate and/or weak magnetic 
field strength, i.e. radiative efficiency is very small. 

In the present study we calculated the smeared im- 
ages by the convolution of the visibility of the theo- 
retical image and the Fourier component of the elliptic 
gaussian beam. When actual observational data are ob- 
tained, however, radio observers usually make images by 
other traditional procedures, such as the method based 
on the CLEAN deconvolution algorithm (see Sec 11.2 of 
TMS and references therein). We can, in principle and 
if successful, produce similar images to what we calcu- 
lated here by the CLEAN deconvolution algorithm, but 
we should be aware in actual data analyses that spuri- 
ous artifacts or some spurious structures are often gen- 
erated in the image during the course of the CLEAN 
deconvolution algorithm (see discussion in TMS). When 
the CLEAN deconvolution algorithm is not successful, 
therefore, the effects of the sidelobe pattern cannot be 
completely removed; i.e. deconvolution errors remain. If 
this is the case, the calculated images do not reproduce 
the true images, and so we cannot correctly extract sci- 



entific information from the data. In order to see this 
explicitly, we show in Fig [7] one extreme case; an image 
created by the dirty beam for model A(b)p05. In this fig- 
ure we can see spurious peaks asymmetric structures in 
the case of the dirty beam. Further we cannot see a cen- 
tral dark region around the black hole in the dirty beam 
case, although it is clearly seen in the case of elliptical 
Gaussian beam. In the analysis of the emission from a 
luminous region around a black hole, it might be difficult 
to make a correct clean beam map, since so far there are 
no good examples reported nor successful observations 
have been made. This is one example of possible diffi- 
culties in the data analysis of the observation of a black 
hole shadow in future. In order to reduce such spurious 
artifacts which might be created in images during the 
process of the CLEAN deconvolution method, other data 
analysis methods have been proposed; e.g. the maximum 
entropy method (MEM, see, Sec 11.3 of TMS). It will be 
interesting and important as future work to investigate 
the validity and the applicability of the data analysis 
methods, such as the CLEAN deconvolution algorithm 
and the MEM, in preparation for future intcrferometric 
observations of the gas flow around black holes. 

5. CONCLUSIONS 

To summarize, we investigate the observational feasi- 
bility of the central structure of the accretion disk, in- 
cluding the black hole shadow, with the space- VLBI ob- 
servations at 22 GHz and 43 GHz by solving the radiative 
transfer equation in the Kerr space-time and adopting 
the accretion flow model proposed for M87. Although in 
this paper we mainly used the parameters of the ASTRO- 
G as an example, the results obtained in this paper can 
be applied to other similar space- VLBI missions. In this 
study, we have obtained following results: 

1 . Our simulations have demonstrated that the asym- 
metry of the intensity map caused by the disk's 
rotation can be detected by the space- VLBI obser- 
vations both at 22 GHz and 43GHz. That is, the 
position of the black hole can be determined and 
this gives the additional evidence for the existence 
of the supermassive black hole in M87. 

2. In order to achieve these scientific goals and those 
described below, the signal-to-noise ratio 7?-sn of 
the observations should be larger than at least 
10 and desirably 50-100 for the observation at 43 
GHz. If the nominal value of the SEFD of 
the ASTRO-G satellite at 43 GHz is achieved, i.e. 
^-satellite ^ 2 8000 Jy, we can expect such SN ratios, 
while for the worst case (i.e. S| atcllitc ~ 190000 
Jy), such SN ratios cannot be expected. At 22 
GHz, the required SEFD to detect the asymmet- 
ric profile is larger than the planned value of the 
VSOP-2/ ASTRO-G by several factors even for the 
worst SEFD value (i.e. S^ tellite ~ 8200 Jy). The 
results at 22 GHz will be applied to the observa- 
tions by the RadioAstron. 

3. From the asymmetric profile, we can put a con- 
straint on the size of the spatial region containing 
the black hole in M87 within 30 GM /c 2 which is 
closer to the event horizon than ever before. If not, 
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the viscosity mechanisms usually proposed by the 
accretion disk theory and the MHD/kinetic sim- 
ulations should be largely altered or be strongly 
constrained. 

4. We have shown that in the cases that the appar- 
ent size of the gravitational radius is r g /D > 15/Lt 
arcseconds (model A or B in Table 1) or that the 
electron temperature profile is not so steep (i.e., 
V ^ 0.7), the under-luminous region around the 
central black hole can also be detected by the space- 
VLBI observations at 43 GHz. Such a relatively flat 
(p 0.7) electron temperature profile is predicted 
by the accretion disk theory and the three dimen- 
sional MHD simulations. If the under-luminous re- 
gion is detected, this gives a very strong evidence 
of the existence of the black hole shadow (or black 
hole silhouette) around the central black hole and 
become the first direct test of the general relativity 
in the strong gravity region around the supermas- 
sive black hole. 



5. We have also shown in the case of r- m = nsco that 
the effects of the black hole spin can also be de- 
tected by the space- VLBI observations at 43 GHz If 
such observation will be perfor med, the spin value 
proposed by the past studies (|Wang et al.l 120081 : 
iLi et al.ll2009D will be independently checked and 
the observation by the space- VLBI observations at 
43 GHz gives more direct evidence of the black hole 
rotation. 
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